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ABSTRACT 

We present a program to study the galaxy luminosity function (GLF) of the Leo I and Coma I Groups 
at ~ 10 Mpc. We have surveyed over seven square degrees in Leo I and ~11 square degrees in Coma I. In 
this paper, we detail the method we have developed and implemented for identifying on morphological 
grounds low-surface-brightness, Mr < — 10 dwarf galaxies at a distance of 10 Mpc. We also describe 
extensive Montecarlo simulations of artificial galaxies which we use to tune our detection algorithms and 
evaluate our detection efficiency and parameter recovery as a function of hr(0) and Rt- We find for a 
sub-set of our Leo I data that at the 90% completeness level we can detect dwarfs comparable to Antlia 
and Sculptor. Finally, we describe preliminary follow-up observations which confirm we are detecting 
dwarf spheroidals in Leo I at 10 Mpc. 

Subject headings: galaxies: dwarf — galaxies: luminosity function, mass function — techniques: image 
processing 



1. INTRODUCTION 

Knowledge of the galaxy luminosity function (GLF) 
is of fundamental importance to extragalactic astronomy 
and cosmology. Studies of the GLF in different density 
environments and at different redshifts provide us with 
the clues to answer questions about how galaxies were 
formed and how their properties and number counts have 
evolved. The local GLF is particularly important for two 
reasons. First, it is virtually impossible to observe faint, 
low-surface-brightness objects at cosmologically interest- 
ing distances. Our understanding of the luminosity density 
of the Universe must therefore come from observations of 
distant, high- luminosity galaxies in combination with an 
extrapolation to fainter objects based on measurements of 
the local GLF (while keeping in mind possible evolutionary 
effects). Second, in order to interpret faint galaxy counts 
in the context of evolution of galaxy properties and to map 
distant populations to local counterparts, the z = GLF 
is necessary for establishing a baseline of current-day prop- 
erties and space densities for galaxies of all luminosities. 

1.1. A Brief History of GLF Measurements 

Because of the importance of the GLF, deriving it in 
different environments has been the subject of many inves- 
tigations in the past two decades. There have been many 
approaches taken to estimate the GLF for the field, in 
the vicinity of luminous galaxies, in groups, in the nearby 
Virgo and Fornax clusters and in more-distant rich clus- 
ters. We make no attempt to give a complete overview 
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of the field here. Pritchet & van den Bergh (1999) have 
a concise summary of recent studies of the GLF in dif- 
ferent environments and Trentham (1998b) reviews GLF 
determinations in somewhat more detail. Impey & Bothun 
(1997) also review recent determinations of the GLF and 
have extensive discussions of issues specifically related to 
selection effects in studies of low-surface-brightness galax- 
ies. 

We use the standard Schechter function (Schechter 
1976) to describe the GLF. The number of galaxies per 
unit volume (<p) at luminosity L ± dL/2 is given by 

(f>(L)dL = (t)*(L/L*) a e- L/L *dL (1) 

where L* is the luminosity of galaxies at the 'break' be- 
tween the exponential form at high L and power law at low 
L and <fi* is the space-density normalization. The charac- 
terization of the GLF to the lowest luminosities — dwarf 
spheroidal (dSph) galaxies like Draco at Mr ~ —9 — is 
likely near complete only for the complement of compan- 
ions to the Milky Way. Yet for the Local Group (LG), M31 
dSph galaxies are still being discovered {e.g. Armandroff, 
Jacoby, & Davies 1999) as are intra-group dSph (e.g. Lav- 
ery & Mighell 1992; Karachentseva & Karachentsev 1998). 
With the currently known members, the Schechter param- 
eters for the LG are M* ~ —21; a ~ —1.1 for the i?-band 
(Pritchet & van den Bergh 1999). 

Beyond the LG, GLF studies have been carried out for 
the field based on redshift catalogues and careful analysis 
of the selection function for the catalogues (e.g. Love- 
day et al. 1992; Marzke et al. 1994; Lin et al. 1996; Zucca 
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et al. 1997; Marzke et al. 1998), in nearby clusters and 
groups (e.g. Sandage, Binggeli, & Tammann 1985; Fer- 
guson & Sandage 1990; Zabludoff & Mulchaey 2000) and 
in more distant clusters (e.g. Bernstein et al. 1995; Tren- 
tham 1998a; de Propris & Pritchet 1998). Despite this 
rich history of measurements, there is no consensus on the 
nature of the faint-end slope of the GLF and how it may 
depend upon environment and morphology. Recent stud- 
ies in clusters have measured a steep faint-end slope (c.f. 
Trentham 1998a; Hradecky et al. 2000), suggesting that 
higher density regions have a higher dwarf-to-giant ratio, 
yet some studies of the low-density field GLF have found 
similarly steep results (c.f. Zucca et al. 1997; Loveday 

1997) . There is also evidence that the faint end varies with 
Hubble-type (e.g. Jerjen & Tammann 1997; Marzke et al. 

1998) , color (e.g. Metcalfe et al. 1998), and emission- 
line strength (e.g. Christlein 2000), and in some cases 
is contrary to what is found in the LG for similar envi- 
ronments. For example, the study of Marzke et al. (1998) 
includes both LG-density groups and the field and reports 
a relatively flat GLF for elliptical and spiral galaxies and 
a steeper GLF for peculiar and irregular galaxies, whereas 
in the LG the irregular GLF is relatively flat. 

The majority of the studies outside the LG to date are 
restricted to galaxies brighter than Mr ~ —16, which only 
begins to probe the faint end of the GLF. Exceptions to 
this are (curiously) some of the studies of relatively dis- 
tant clusters where the claimed low-luminosity limits are 
Mr ~ -12 (e.g. Coma, Bernstein et al. 1995). These 
studies have generally taken a statistical approach in sub- 
tracting off background and foreground contamination and 
have pushed the faint limits of their GLFs to the limits of 
the photometry. Potential problems with this approach 
(particularly at faint levels where the surface density of 
distant galaxies is high) are projection effects (Valotto, 
Moore, & Lambas 2000) and the existence of large-scale 
structure giving rise to significant variations in the back- 
ground counts as a function of position. Most studies con- 
centrating on nearby groups and clusters (necessarily ex- 
tending over large areas on the sky) have identified low- 
surface-brightness (LSB) dwarfs on photographic plates, 
usually followed up with CCD observations (e.g. in Virgo 
and M81; Impey, Bothun, & Malin 1988, Caldwell et al. 
1998). Only recently have wide-field CCDs been available 
to cover similarly large expanses of the sky (c.f. Carrasco 
et al. 2001). 

1.2. A New Survey 

We have undertaken a program to measure the GLF to 
very faint levels in nearby (~ 10 Mpc < d < 15 Mpc) 
galaxy groups using i?-band MOSAIC CCD images for 
the initial survey and a number of different techniques to 
establish group membership of candidates on a galaxy-by- 
galaxy basis. Our primary goal is to measure the z = 
GLF to luminosity and completeness limits similar to that 
achieved in the LG. This will allow a test of the variation 
of the GLF between groups, decrease the shot noise in the 
GLF by combining samples or simply having a larger sam- 
ple of galaxies than is present in the LG. Furthermore, it 
will provide a large sample of low-luminosity galaxies that 
are near enough to be amenable to fairly detailed studies 
of individual galaxy properties. 



At these distances, dwarf galaxies, at least those repre- 
sented in the LG, have apparent sizes > 10" at our surface- 
brightness limit (50% completeness) of fiR ~ 26 mag/// 2 . 
It may therefore be possible to identify dwarf-galaxy mem- 
bers purely on morphological grounds. The imaging data 
of our survey were taken with the MOSAIC camera on 
the KPNO 0.9-meter telescope (described more fully in §2) 
with the engineering-grade CCDs. Due to the poor blue 
response of these original chips, we imaged our survey with 
the Cousins R filter. Figure 1 shows the positions of LG 
dwarfs as they would be seen at 10 Mpc in a plot of total 
magnitude vs. central surface brightness. LG dwarf prop- 
erties are taken from Mateo (1998) and references therein, 
adopting (B — R) ~ 1.3. The dashed lines show apparent 
sizes of galaxies with exponential profiles in this plane at a 
limiting surface brightness isophote of [ir = 26.7 mag/// 2 . 
Most of the LG galaxies have diameters at this isophote 
level that are greater than 10". 

We describe in detail our procedures for identifying in- 
trinsically faint dwarf galaxies in the following sections. 
We note here that there are not very many types of ob- 
jects known with the low central surface brightness level of 
a typical dSph. Distant high-luminosity galaxies are sub- 
ject to (1 + z) 4 surface brightness dimming and enter into 
our regime of interest, but these objects will typically be 
apparently small — a few arcseconds or less at our isopho- 
tal limit. Background large low-surface-brightness galaxies 
are potential sources of confusion with our method (Dal- 
canton et al. 1997) as are very distant galaxy clusters for 
which the main signature is a large, low-surface brightness 
enhancement due to the overlapping envelopes of the clus- 
ter galaxies (Zaritsky et al. 1997). However, with a some- 
what conservative isophotal diameter limit of 10" these 
objects should not be a significant contaminant. We will 
carry out group-membership studies (see §5.2) for at least 
a subset of our morphologically-identified candidate mem- 
bers. 

The initial groups in this survey are the Leo I group 
(NBGC 15-01 or LGG 217, containing NGC 3379 and 
NGC 3368; Tully 1988; Garcia 1993) and the Coma I 
'cloud' (NBGC 14-01 or LGG 279, containing NGC 4278 
and NGC 4314; Tully 1988; Garcia 1993). In this paper we 
describe our algorithms for searching the frames for candi- 
date group members, techniques for deriving photometric 
and structural properties of the candidates and simulations 
used to tune the search technique parameters and evalu- 
ate completeness and measuring errors. We also present 
examples of the initial findings. In a subsequent paper we 
will present the details of the observations, our candidate 
lists, charts, and the GLF. 

2. SURVEY IMAGING 

This program begins with direct imaging using the MO- 
SAIC Camera on the 0.9-meter telescope at KPNO. This 
camera gives a field of view 1° on a side and is comprised 
of eight 2048 x 4096 CCDs, with sampling (unbinned) of 
0.42 arcsec/pixel. The details of the observations for each 
group will be presented in the subsequent data papers. 
In short, nine fields were observed in the Leo Group and 
twelve fields in the Coma I Group. For most of the Leo 
Group observations the original engineering-grade CCDs 
of the MOSAIC were used. For each field, five 900-sccond 
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exposures were made through the Cousins R filter using 
the standard dithering routine with telescope offsets of 
~ 40" between the images. There were significant gaps 
between the CCDs (~ 21" N-S, - 15" E-W ) in the MO- 
SAIC and numerous cosmetic flaws in the engineering- 
grade CCDs. The purpose of the in-field dithering was 
to allow removal of these artifacts in the process of com- 
bining the five frames into a single image. In retrospect, 
we should have had a larger number of dithered point- 
ings to better remove artifacts and improve our ability to 
flat-field on intermediate scales. 

Our procedures for bias and dark subtraction, flat- 
fielding and combining the frames followed the steps de- 
tailed in the Guide to the MOSAIC Data Handling System 
(Valdes 1998). To carry out this initial image processing 
we used the MSCRED package within IRAF 3 . The prelimi- 
nary processing of MOSAIC frames is somewhat compli- 
cated and will be discussed in detail in a subsequent paper. 
To do photometric calibration of the stacked MOSAIC im- 
ages, we independently observed our Leo I and Coma I 
fields in B and R with the Lick Observatory 1-m Nickel 
telescope. The calibration is good to ±0.04 mag in R, but 
when we include the uncertainty in the color-term from 
our estimates of the B — R colors of the detected objects, 
the final galaxy photometric error is ±0.07 mag. Figure 2 
highlights one of the nine fields imaged in the Leo I Group, 
with the 0.9-m fields overlaid upon a Digital Sky Survey 
image 4 . The ~ 3°x ^3° imaging area encompasses the 
1?1 virial radius of the group (Tully 1987). We analyze 
the indicated field (Field 7) to demonstrate our search and 
simulation algorithms. Our coverage of the Coma I group, 
while larger in extent (~ 11 deg 2 ) covers a smaller fraction 
of the extensive group (> 25 deg 2 ). 

3. THE ANALYSIS PIPELINE: DETECTING FAINT 
GALAXIES 

We have developed a three-step process for detecting 
dwarf galaxies at 10 Mpc. With the first, we detect high- 
surface-brightness (HSB) galaxies through a standard k-a 
thresholding technique using the SExtractor object detec- 
tion and classification software (Bertin & Arnouts 1996). 
In the second step, we use a search method we have opti- 
mized to find LG-like dwarf galaxies that would typically 
fall below the detection limits of most standard meth- 
ods. We use a matched-filter technique for which the HSB 
galaxies and bright stars are masked, and the resulting 
image is convolved with an exponential filter representa- 
tive of the structure of a LG dwarf at a distance of 10 
Mpc. This second step of our detection process is based in 
part on the technique developed by Dalcanton (1995) to 
detect large low-surface-brightness (LSB) galaxies in the 
field from drift-scan data. In the third step of our analysis 
pipeline we characterize the detections and create a reli- 
able candidate list from which to do further analysis and 
observations. In this section, we describe our detection 
methods and illustrate them in Figures 3-6. 



3.1. Step 1: Traditional Detection Method 

Before object detection, each image was pre-processed 
to remove the brightest objects from the image. We mod- 
eled the largest, most luminous galaxies in the field using 
the IRAF task ELLIPSE, and subtracted off their elliptical 
light components. This revealed objects near the line of 
sight of the large galaxies, as well as remaining structure 
in the galaxies. We then removed the outer light from all 
stars brighter than m# ~ 10 — a somewhat arbitrary cut 
that removed the largest stars with diffraction spikes, but 
didn't present an unreasonable number of stars to fit by 
hand. The saturated centers and diffraction spikes of the 
stars were masked with an IRAF routine that flags objects 
brighter than a specified threshold, and these masked ob- 
jects were then replaced with a median-smoothed approx- 
imation of the background from SExtractor plus artificial 
noise added with the IRAF task MKN0ISE. We modeled 
the resulting smoothed star image with ELLIPSE, and sub- 
tracted the modeled light from the original, unmasked im- 
age. 

Once we rid our image of diffuse light from bright ob- 
jects, we ran SExtractor in 'single image mode', where 
detection, photometry, and shape measurements were all 
done on the same input image. We set a 70-pixel minimum 
contiguous area for detection and a detection threshold of 
1.5-a s k y above the sky. We required SExtractor to create 
noise maps of the image internally, and weight detections 
according to the noise characteristics of the section of the 
image where the detection was made. For star /galaxy 
separation, we used the object classifier of SExtractor, 
CLASS_STAR, and kept only objects with CLASS_STAR < 0.1 
and semi- major axis of the SExtractor ellipse fit, A_W0RLD 
> 1", which corresponds roughly to an isophotal radius 
of 3" (we found that for the smallest objects, SExtrac- 
tor could no longer distinguish stars from galaxies). With 
the remaining objects, we fit radial profiles with ELLIPSE 
in IRAF, inputting the measured SExtractor shape pa- 
rameters as initial fit values. We fit the one-dimensional 
isophotal profile of each galaxy output from ELLIPSE, as- 
suming an exponential profile. We then characterize our 
detections from the traditional method by adopting the 
central surface brightness as extrapolated from the radial 
profile fit, and the magnitude as measured by SExtractor, 
MAG_BEST. Of these objects we retain those with a limit- 
ing isophotal size at im m = 26.7 of @u m > 10" assuming 
an exponential profile with zero ellipticity (Allen & Shu 
1979). Further object selection was done through morpho- 
logical membership selection on a galaxy-by-galaxy basis, 
examining each object by eye (described in §5). We used 
simulations (see §4) to tunc all detection parameters used 
in this step. 

To demonstrate the first step of our galaxy detection 
process, we show in Figure 3 a 15' x 15' sub-section of our 
Field 7 image, on which we have placed simulated images 
of LG galaxies as they would appear at a distance of 10 
Mpc (the galaxy simulation procedure is described in §4). 



3 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science Foundation. 

4 Based on photographic data obtained using Oschin Schmidt Telescope on Palomar Mountain. The Palomar Observatory Sky Survey was 
funded by the National Geographic Society. The Oschin Schmidt Telescope is operated by the California Institute of Technology and Palomar 
Observatory. The plates were processed into the present compressed digital format with their permission. The Digitized Sky Survey was 
produced at the Space Telescope Science Institute (ST Scl) under U. S. Government grant NAG W-2166. 
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We have modeled all dSphs as exponentials. Note that 
at this distance Leo II, Tucana, And V, and Carina are 
extremely difficult to detect by eye. In Figure 4, we show 
the results of Step 1, where an aperture is drawn around 
every object detected that is also classified as a galaxy by 
SExtractor (CLASS_STAR < O.f, A.WORLD >1"). Many of 
the artificial LG dwarfs are detected; however, the lowest- 
surface-brightness objects, And II, Tucana, And V and 
Carina, axe fainter than the 1.5-<r s fcy 

detection threshold. 

3.2. Step 2: Optimized Detection Method 

With the second method, we essentially removed all ob- 
jects detected by the first method, and convolved the re- 
maining background with a smoothing kernel to bring out 
LSB features. To remove all bright objects, we masked 
all HSB objects above a 2-a s k y surface-brightness thresh- 
old. We replaced the masked pixels in our field with a 
single background value plus noise added with MKNDISE 
(appropriate to the field). We then convolved the masked 
image, using the IRAF task CONVOLVE, with a kernel op- 
timized to enhance fluctuations of the size and shape of 
LG-sizcd dwarf galaxies at 10 Mpc. In this example we 
used an exponential kernel with a 5" scale length, as dSphs 
are known to be reasonably well-fit by exponential profiles 
(7 = 7 e- (r/a) ; c.f. Faber & Lin 1983; Vader & Chaboyer 
1994), and 5" is a characteristic exponential scale length 
for LG dSphs if moved to 10 Mpc (Mateo 1998). We find, 
however, that our detection rate is not overly sensitive to 
the shape of the convolving kernel. 

We ran SExtractor in 'double image mode' on the result- 
ing convolved images, where the detection and shape mea- 
surements were performed on the convolved image, but 
the corresponding photometric measurements were done 
on the original image. We set the detection limit to re- 
quire a minimum area of 100 pixels above a 5.5-er s fc y sur- 
face brightness threshold. The simulations in §4 were used 
to optimize these parameters. SExtractor's photometry in 
double- image mode, however, tends to underestimate mag- 
nitudes and sometimes central surface brightnesses due 
to its inability to de-blend neighbor objects on the sec- 
ond (measurement) image. De-blending occurs here on 
the convolved image, where bright neighbors have been 
masked and so not included in de-blending. For final pa- 
rameter measurement we fit radial profiles for each opti- 
mized detection, as discussed below. 

There are a small number of spurious detections from 
this step, as the convolution process not only enhances 
features of a certain scale but also smoothes them in 
such a way that SExtractor's neural network training for 
CLASSJ3TAR can be fooled. However, most brightness en- 
hancements detected in the convolved image can be culled 
by hand, as most of these false detections are obviously 
poorly masked stars and brighter galaxies. An example 
of such detections can be seen in Figure 5, where we show 
the masked, convolved image that corresponds to Figure 3. 
Here, all SExtractor detections from Step 2 are indicated: 
the obvious masking remnants detected are marked with 
crosses, and the candidate detections are marked with el- 
liptical apertures. 

In this example, all artificial LG dwarfs are detected, ex- 
cept Leo II. Note the large size of the elliptical apertures, 
which are roughly equivalent to the isophotal radii, for the 



input galaxy detections. We use the large apparent size of 
such dwarf galaxies to discriminate against more distant 
and typically apparently smaller galaxies. 

One challenge for this search method is that the masking 
process may mask the higher-surface-brightness centers of 
the sort of very faint and extended galaxies we hope to 
detect. The convolution process may then smooth the re- 
maining light, but the resulting convolved object will have 
a fainter peak luminosity and perhaps fall below the de- 
tection threshold. An example of this in Figure 5 is Leo 
II, which was detected in Figure 4, but is masked and 
undetected by the optimized method. One possible prob- 
lem would be the case where an object falls just below the 
threshold of the traditional method, but is also masked and 
overlooked by the optimized method. To guard against 
these losses, we set the masking threshold (2-osky) to be 
somewhat higher than the detection threshold for the tra- 
ditional method (1.5-a s k y ), so that even though some ob- 
jects will be partially masked in the optimized method, 
they are likely to have already been detected in Step 1. 
In this way, we have a number of intermediate-surface- 
brightness objects that are detected twice, but we mini- 
mize the chances of complete non-detections of such ob- 
jects in this middle ground. 

Objects detected via both methods, or by the opti- 
mized method alone, have an extended LSB component 
that falls below the traditional surface brightness thresh- 
old. For these objects, the traditional detection aperture 
misses this extended light, underestimating the total flux. 
Furthermore, the SExtractor magnitude measured by the 
optimized method tends to be skewed by bright neigh- 
bors that may fall within the measurement aperture, as 
mentioned earlier. We avoid these problems by measuring 
the isophotal radial profile with ELLIPSE, again fitting the 
profile with an exponential. We then characterize these 
detections with the fitted central surface brightness, 
as in the traditional method, but for the magnitude, we 
adopt the total magnitude calculated analytically from /!q* 
and the fitted exponential scale length, a fit . 

In Figure 6, we plot the example detections from Fig- 
ures 4 and 5 on the Rt vs. /iq plane. Here and in the 
following figures, Rt is the total magnitude as measured 
by each detection method. Detections from the traditional 
method are indicated as crosses, where we plot the total 
magnitude measured by SExtractor (MAG_BEST) and the 
central surface brightness parameters from the fitted pro- 
file. As in seen in §4.3, for these higher-surface-brightness 
galaxies, the profile fit does a good job estimating the cen- 
tral surface brightness, but uncertainties in the slope of 
the profile yield larger uncertainties in the calculated total 
magnitude than we see in the SExtractor magnitude mea- 
surement. We plot optimized detections as stars, where 
Rt is now the total magnitude calculated analytically from 
/Ltg* and the fitted exponential scale length, a fit . The ar- 
tificial input galaxies that were successfully detected are 
circled. Even though all the included objects are detected, 
Carina and And V represent the very limits of our de- 
tection method as will be seen in the following section. 
As mentioned previously, the proximity of the Leo I and 
Coma I groups makes dwarf galaxies large in angular size. 
This can be seen in Figure 6 in the separation of smaller, 
more compact background objects (the dense concentra- 
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tion of crosses centered around Rt ~ 19.5, /zo ~ 22) and 
LG galaxy properties. Membership characterization will 
be discussed more fully in §5. The input galaxies (circled 
objects) all are measured with parameters comparable to 
their input values and so fall along the main locus of LG- 
like galaxies. Our ability to measure these parameters is 
explored with the simulations, as described in the next 
section. 

4. SIMULATIONS 

Two factors affect our ability to measure the GLF: 
(1) contamination from non-group members and artifacts, 
such as noise peaks or flat-fielding errors mistakenly se- 
lected through the optimized detection method, and (2) 
incompleteness in our detection method. We address the 
first on a galaxy-by-galaxy basis, guided by measured 
properties, morphology, and follow-up observations. The 
second we address via Montecarlo simulations in which 
we added artificial galaxies to the data frames and at- 
tempted to recover them with our detection processes de- 
scribed in §3. The final luminosity function will be 
membership-corrected for each of these effects. The final 
galaxy counts are the observed counts of member galaxies, 
n bs — n b s (RT , Ho) , weighted by the detection efficiency 
as determined from the simulations, /(Rt, Ho), where Rt 
and no are the input total magnitude and central surface 
brightness, respectively, of the artificial galaxies. Thus, 
the final galaxy counts for the GLF will be N(Rt,Ho) = 
nobs/ f ■ We present a subset of these simulations here; the 
full simulations results will be presented in a subsequent 
paper. 

4.1. Simulation Procedure 

In order to understand the complex selection function 
and to tune our detection parameters for both steps in our 
detection method we have carried out extensive artificial- 
galaxy simulations. We create artificial galaxies covering 
the full range of the LG-dwarf magnitudes and central sur- 
face brightnesses (as seen at 10 Mpc), degrade them by the 
seeing and noise properties of the real data, add them to 
the real data frames and then detect them by the method 
described in §3. We calculate the detection efficiency as 
/(Rt, Ho) = n rec /n a dd, where n rec is the number of galax- 
ies recovered and n a dd is the number of galaxies added for 
a given bin. The artificial galaxies populate the Rt vs. ho 
plane spanning the ranges 14 < Rt < 22, 22 < Ho < 26, 
roughly corresponding to the range wherein we would see 
LG-counterpart dSphs. 

The artificial galaxies we created with the IRAF task 
MKDBJECT in the same manner as those shown in Fig- 
ure 3. With MKDBJECT we input the magnitude zeropoint 
(which includes the color term for an adopted dwarf color 
of B — R = 1.3), read-noise, and seeing of our observa- 
tions, and we specify an exponential radial profile, total 
magnitude, characteristic scale length, axial ratio, and po- 
sition angle. The routine approximates the profile analyti- 
cally and artificially truncates the profile at a cut-off radius 
where the profile intensity is 10~ 5 times the peak inten- 
sity (~ 12.5 mag below the peak). The dwarf galaxies we 
simulate are not significantly affected by this truncated dy- 
namic range. Again, we adopt a purely exponential profile 
as it is known to be a good model for faint dSph galaxies 



and we include a uniform range of axial ratios to model 
the galaxies' intrinsic ellipticities. LG compact dwarf el- 
lipticals like M32 follow a de Vaucoulcurs i? 1 / 4 -law, but 
their high central surface brightnesses (as can be seen by 
M32's position in Figure 1) will make them easy to detect 
at these distances. Thus, we do not include these profile 
types in our simulations. There have long been discussions 
in the literature comparing the three-parameter King pro- 
file with the two-parameter exponential profile as the most 
appropriate for early-type dwarfs (see Ferguson & Binggeli 
1994 and references therein), with exponentials being very 
similar to a King profile with large core radius and low con- 
centration. Fainter dwarfs (Mr < —17) appear to be well- 
fit by an exponential at all radii yet brighter (Mr < — 17) 
dwarfs, for which we see some profile flattening (e.g. For- 
nax dwarfs; Bothun, Impcy, & Malin 1991) or excesses in 
the inner regions (Caldwell & Bothun 1987), benefit from 
the more flexible King profile. We adopt an exponential 
profile for our fake galaxies for several reasons. First, a 
two-component model simplifies the parameter space we 
explore with our simulations. Second, deviations between 
these profiles in the core regions are minimized by the see- 
ing convolution of our simulations at 10 Mpc. Finally, 
our simulations are most critical for the lowest luminos- 
ity regime (Mr > —14), where the exponential profile is 
found to be the most prevalent. 

For these simulations we neglected any internal struc- 
ture, patchincss due to dust or HII regions and thus inter- 
nal absorption. For dSph galaxies, this is appropriate. To 
a first approximation, this also tells us something about 
our i?-band detection efficiency for dwarf irregulars, which 
typically have smoother and more extended disk of old, red 
stars underlying a younger and bluer stellar population 
that is well-fit by an exponential (Irwin & Hatzidimitriou 
1995). However, in order to fully probe the effects of sub- 
structure on detectability, we will add high signal-to-noise 
images of real LG dwarf irregulars (e.g. GR 8, IC 10) to 
our frames, scaled to a distance of 10 Mpc with compa- 
rable noise and seeing to our data, and will attempt to 
detect them with our method. This will be included in a 
subsequent paper. 

We have randomly generated galaxy parameters uni- 
formly in bins of size 0.5 mag/// 2 in ho and 0.5 magnitudes 
in Rt, according to the appropriate radial profile. We then 
randomly distributed each galaxy in the image frame with 
a spatial density that allowed for a possible < 1% overlap 
of the area in input objects (assuming a maximum isopho- 
tal size at the 27 mag/// 2 isophote). Preliminary simu- 
lations were first performed by adding exponential-profile 
galaxies with 17 < R T < 22, 176 < a < 6'.' 3 to a purely 
artificial noise frame, mimicking the noise characteristics 
of the CCD, in order to tune the simulation and detec- 
tion algorithms. For the traditional method, we tuned the 
detection parameters solely to maximize the probability 
of detection. We found that a detection threshold of 1.5- 
o- s k y vs. 2-a s ky made little difference, and so we chose the 
more generous threshold. We reduced the minimum de- 
tection area (DETECT_MINAREA) based on the size of artifi- 
cial LG galaxies moved to 10 Mpc, adopting the minimum 
isophotal area above the detection threshold of the faintest 
LG-likc galaxies we expect to detect with the traditional 
method. Lastly, these simulations also helped us to tunc 
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the background mesh sizes (BACK_SIZE), which governs 
how large a background area is effectively smoothed. We 
found the size did not significantly change the detection 
efficiencies. 

The optimized detection parameters were also tuned 
from simulations with these purely artificial noise frames, 
which we found to be unrealistically flattened images com- 
pared with simulations done on the real data images. How- 
ever, the general detection parameters could still be opti- 
mized with these more idealized simulations. We found 
that extremely low detection thresholds identified many 
random noise peaks in the convolved image as real detec- 
tions. These preliminary simulations also tested the mini- 
mum area parameter, DETECT_MINAREA, and the kernel size 
of the convolution filter. With a large DETECT J4INAREA of 
300 pixels we detected objects only along a narrow central 
surface brightness regime, and lost a significant number of 
small objects. Using a DETECT_MINAREA of 100 pixels, we 
recovered almost all input objects in the parameter space 
tested. Changing the kernel size from a 5" exponential to 
a 10" exponential did not change the results dramatically, 
except insofar as the smaller kernel pushed the limiting 
central surface brightness a little fainter; thus, we adopted 
the smaller filter size. 

Ultimately, we expanded the simulations, adding the ob- 
jects to the real data frame for Field 7. The image has 
had the brightest galaxies modeled and subtracted, and 
the halos of light surrounding the brightest stars modeled 
and subtracted according to the procedure described in 
§3.2. Distributing galaxies in the real image probes not 
just the efficiency of the detection parameters, but also 
loss of chip area due to very bright objects and detection 
difficulties due to background inhomogeneities. We find 
that the latter difficulties with the background dominate 
the detection process in the optimized method. 

With real data images we further tuned our optimized 
search method to avoid possible false detections. False 
detections can arise in the masked and convolved image 
from peaks in the noise and from flat-fielding artifacts 
(e.g. edges from the individual MOSAIC CCDs). Follow- 
up imaging of candidates and use of detectors with bet- 
ter cosmetics than the MOSAIC engineering-grade chips 
will ultimately eliminate almost all of these false detec- 
tions. In the meantime, we can characterize the expected 
contamination from false detections, and tune our detec- 
tion threshold to exclude them, using a negative version 
the current MOSAIC image. We first inverted the masked 
and background-subtracted version of the image, and then 
added the background back in. The resulting image is then 
convolved with the same kernel, and the optimized detec- 
tion method applied. By inverting the image, we insure 
that the only possible detections are artifacts and not real 
objects. False detections from noise should happen equally 
often for the positive as for the negative image. Similarly, 
the MOSAIC CCD edge artifacts are the result of very 
slight illumination mismatches — both positive and nega- 
tive — between the dithered images, so they create both 
positive and negative sharp edges in the convolved image. 
Other, more bizarre flat-fielding errors create structure 
that is also typically symmetric around the background 
level. Thus, on average, the inverted image should result 
in the same number of false detections from these types of 



contamination sources. For Field 7, we set our detection 
threshold for the optimized method to b.5-o~ s k y , at which 
level we find zero possible false detections that cannot eas- 
ily be identified with an obvious flaw. 

The total number of objects simulated in a parameter 
bin, n a dd(RT, Mo), is determined by the criterion that the 
error in the final counts N(Rt,Ho) contributed from the 
uncertainty in / is small compared to the error contributed 
by Poisson variance in the observed counts, n b s (RT, Mo)- 
We adopt the error model for deriving stellar luminosity 
functions corrected by artificial star experiments of Bolte 
(1989). The uncertainty in a bin is given by 

n bs . (l-j>L 
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We add a preliminary minimum number of 100 galaxies 
for every bin and then run the experiments adding galaxies 
from appropriate bins until Equation 4 is satisfied. 

4.2. Selection Function 

The explored parameter space presented in this paper, 
17 < R T < 22, 22 < u < 27, is a subset of the final sim- 
ulation results to be presented in a subsequent paper, but 
encompasses the transition region between our traditional 
vs. optimized methods for the lowest luminosity dwarfs. 
The completeness fraction from these simulations is plot- 
ted in greyscale in Figure 7 as a function of Rt and Uq, 
where the greyscale levels correspond to 10% (lightest), 
30%, 50%, 70%, and 90% (darkest) completeness. The 
first two panels, Figures 7a and 7b, show / as determined 
independently by the two detection methods, traditional 
and optimized, respectively. The fact that the optimized 
method complements the traditional method can be seen 
in these contours, where the optimized method pushes the 
completeness envelope to fainter surface brightness limits. 
The composite completeness factor combining both meth- 
ods is represented in Figure 7c. LG galaxies if seen at a 
distance of 10 Mpc are plotted for comparison as filled cir- 
cles, showing that at the ^90% completeness level we can 
detect dwarfs comparable to Antlia and Sculptor and at 
the ~50% completeness level we can detect dwarfs simi- 
lar to Tucana and Leo II. The complexity of the selection 
function shows that the completeness of our survey can- 
not be characterized by just a magnitude limit or a sur- 
face brightness limit, but requires combining both effects. 
We can say that the faintest central surface brightness we 
reach at the > 50% level is around 25.7 mag/// 2 , and the 
faintest magnitude is around Rt — 20. In measuring pro- 
files of our real and artificial galaxies, we robustly detect 
our objects to a limiting surface brightness isophote of 26.7 
mag/// 2 . 

4.3. Parameter Recovery 

The simulated galaxy experiments can also give a mea- 
sure of how robustly we measure structural parame- 
ters. For each object recovered, the radial profile is non- 
intcractively measured during the simulation using the 
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IRAF task ELLIPSE. Wc fit an exponential profile and de- 
termine the extrapolated central surface brightness, hq, 
and the fitted exponential scale length, a. We do not test 
our ability to discriminate between different profile shapes 
here. 

In Figures 8-11, we demonstrate our ability to recover 
the input parameters central surface brightness, total mag- 
nitude, and exponential scale length. Our parameter re- 
covery accuracy depends strongly upon the position of the 
galaxy in Rt vs. hq parameter space, and as expected, 
significantly decreases as we approach our detection limits 
in either axis. Figure 8 shows the difference between input 
and output ^ plotted as a function of input total magni- 
tude, with the left panels showing all the artificial galaxies 
recovered by the traditional method, and the right pan- 
els those recovered only by the optimized method. In the 
vertical direction, we divide the artificial galaxies by input 
central surface brightness: 22 < ^ < 23, 23 < ^ < 24 
and 24 < (i™ < 25 from top to bottom. We can see our 
accuracy in parameter recovery decrease as the detection 
efficiency starts to fall, near Rt — 20. The limits on n a dd 
imposed by Equation 4 can be seen, for example, in pan- 
els a, b, and e, where there is a dramatic increase in the 
number of simulated galaxies for certain bins. The com- 
plementary nature of the two detection methods is also 
implicit in this figure and can be compared to Figure 7. 
Panel a shows the traditional method peaking in efficiency 
for /if," ~ 22.5, 17 < Rt < 19.5, and panel f shows the op- 
timized method peaking for /i n ~ 24.5, Rt — 17.5. We see 
in the right-hand panels that the lowest-surface-brightness 
galaxies detected only through the optimized method have 
a large scatter in A^ - These objects are more difficult to 
fit due to their very low S/N which follows from their very 
low surface brightnesses, especially as we approach our de- 
tection limits {e.g. panel e). 

In Figures 9 & 10, we show the magnitude recovery from 
the simulations, comparing the input magnitude to the 
SExtractor-measured magnitude in Figure 9, and to the 
total magnitude calculated from the fitted parameters in 
Figure 10. In each figure we show the difference in input 
and output total magnitude, ARt = R T ut — Rt > as a func- 
tion of R™, with the same panel divisions as the previous 
figure. First, we recover the expected offset in SExtractor 
MAG_BEST in the left-hand panels of Figure 9. This offset 
results from SExtractor's use of an aperture correction to 
approximate the total magnitude which is known to sys- 
tematically miss some light at the ~ 0.06 magnitude level 
(Bertin & Arnouts 1996). Second, we see in the right- 
hand panels of Figure 9, especially in panel d, evidence for 
the inability of SExtractor to de-blend objects properly 
when run in double-image mode, as discussed previously. 
We see the faintest galaxies tend to be measured brighter 
than they were input, as bright, neighboring objects that 
fall within the photometric aperture are mistakenly being 
included in the measurement. This would also introduce 
larger deviations from the true magnitude for fainter ob- 
jects, and so we see the scatter increase in panels e and f. 
In comparing the results from the traditional method in 
Figures 9 and 10, we see that SExtractor MAGJ3EST does 
a slightly better job of measuring the magnitude than the 
total magnitude calculated from the profile fits, which is 
why we adopt MAGJ3EST for these detections. The larger 



scatter in ARt seen in the left-hand panels of Figure 10 are 
primarily due to difficulties in measuring the exponential 
scale length in the radial profiles, as can be seen in Figure 
11. The considerable scatter in the right-hand panels of 
Figure 10 shows that, even inputting artificial exponential 
profiles, the noise at these low surface brightnesses makes 
automated profile fitting a difficult endeavor. 

In Figure 11 we demonstrate our ability to recover the 
exponential scale length (a) of the simulated galaxies with 
our profile fitting. Here we show the difference in input 
vs. output exponential scale length, Aa = a out — a m , 
as a function of R™. Panels a - f are divided by central 
surface brightness in the same way as the previous fig- 
ure. While the extrapolated central surface brightnesses 
for these galaxies are well recovered, as seen in Figure 8, 
we find that the slopes of the profile fits, which yield the 
exponential scale lengths, are much less well determined. 
An unsuccessful profile fit tends to approach the sky level, 
as background subtraction can be difficult for a low S/N 
object. On average this yields flatter fits for more difficult 
objects, and thus larger a out . We see this in our sim- 
ulations, where the scatter in Aa always tends towards 
a out > a m . Furthermore, the errors in a ont seem to domi- 
nate the calculated total magnitude at all magnitudes and 
central surface brightnesses, as the scatter in Figure 11 to 
preferentially larger a out roughly corresponds to scatter in 
Figure 10 to preferentially brighter R T ut (i.e. ARt < 0). 
Finally, we see again that for the majority of the objects 
detected only through the optimized method (which are 
typically of the lowest surface brightness or lowest S/N), 
a is very difficult to determine. 

Some of our difficulties in measuring the profiles of the 
faintest and lowest-surface-brightness galaxies are miti- 
gated for our real candidate detections by interactive pro- 
file fitting and deeper follow-up imaging in the future. We 
note that most photometric studies of galaxies with sim- 
ilar properties as our artificial galaxies are likely affected 
to similar levels. 

5. PRELIMINARY RESULTS 

To demonstrate the success of our method, we present 
a preliminary sample of objects detected in Leo I. First, 
we show the measured properties of objects found over the 
one-square-degree image of Field 7, and discuss how these 
properties can illuminate these objects' group member- 
ship. Second, we examine several objects from the group 
in more detail, presenting examples of the types of supple- 
mentary observations we will also use to constrain mem- 
bership and contamination. Among these objects wc in- 
clude a newly discovered dSph in Leo I which has compara- 
ble properties to a dSph like Sculptor (Mr = —11.11, /io = 
23.39 R mag/// 2 ). We also show exciting examples of how 
we may use the classic test of resolving our dwarf galaxies 
into stars to establish their membership as done by Baade 
(1944) to confirm M32 and NGC 205 as companions to 
M31, as well as using surface brightness fluctuations (SBF) 
for relative distance determination. Lastly, we show a pre- 
view of the contribution these detections make to the faint 
end of the GLF. 

5.1. Field 7 Objects 
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As mentioned previously, our object search process and 
membership determination are guided by our knowledge 
of objects found in the LG. We know that most dwarfs 
at a distance of 10 Mpc will subtend a large angle on the 
sky (> 20"), plus appear distinctly diffuse (LSB) with a 
roughly flat radial profile. This was shown visually in Fig- 
ure 3. In addition, the more luminous late-type galaxies 
at 10 Mpc will have resolved structure such as spiral arms, 
patchiness and HII regions. Early-type galaxies, the LG 
example being M32, have extremely bright central surface 
brightnesses and steep radial profiles which fall off quickly 
in angular size with increasing distance. These morpho- 
logical criteria, covering the most common morphological 
types we find locally, allow for elementary membership as- 
signment which distinguishes group members from very 
distant background galaxies (c.f. Sandage & Binggeli 
1984). Although these criteria have been very success- 
fuly employed, for instance, for the Virgo Cluster Cata- 
log (Drinkwater et al. 1996), this morphological classifica- 
tion can begin to break down for more nearby background 
galaxies, as well as for more unusual types of galaxies not 
found in the LG. Such types of galaxies are dwarf SO galax- 
ies (Sandage & Binggeli 1984), BCGs (Thuan 1983) and 
their subset, the so-called HII galaxies (Telles & Terlevich 
1997), and large, very low-surface-brightness disks (Impey 
et al. 1988). In these less certain cases, we can begin to 
discriminate non-group-member galaxies by further char- 
acterizing our sample using their measured properties such 
as total magnitude, central surface brightness, and radial 
profile shape. 

We present the properties of our detections from the en- 
tire one-square-degree image of Field 7 in Figure 12, plot- 
ted in the same way as Figure 6. The detections in this 
figure include contamination from background non-group- 
members, which appear most prominently near Rt ~ 
19.5, jLto ~ 22. The crosses show all detections from the tra- 
ditional method after star/galaxy separation (as in Figure 
6). The plotted parameters for these objects (the crosses) 
are the total magnitude as measured by SExtractor and 
the central surface brightness as extrapolated from the ra- 
dial profile fit. Objects detected through the optimized 
method are shown as stars. The plotted parameters for 
these detections are both calculated from the profile fit. 
Objects detected by both methods are connected by a line. 
For comparison we show again LG dwarfs scaled to a dis- 
tance of 10 Mpc, which are marked by filled circles. 

Again we see the effectiveness of this plot as the smaller, 
more compact objects that are mostly distant background 
galaxies separate out from the expected locus of LG-like 
dwarfs. This also follows from the radial profiles, and thus 
morphologies, we expect from known galaxy types, which 
can be qualitatively seen in Figure 1. For example, early- 
type, de-Vaucouleurs-profile galaxies do not populate this 
LSB regime of parameter space shown in Figure 12 un- 
less they are very distant, cosmologically dimmed objects, 
and even then their small angular sizes would make them 
virtually undetectable within our limits. 

The galaxies circled in Figure 12 are described in more 
detail in Figures 13 - 16, where we show a 2'x2' postage 
stamp of the MOSAIC i?-band image for each example 
and its radial profile measured with ELLIPSE. In Figure 
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13a we show one of the larger, higher-SB galaxies in Field 
7. This is a previously catalogued object, CGCG 065-091 
(Zwicky et al. 1968), which has a measured velocity listed 
in NED 5 of 7922 km s" 1 (Fouque et al. 1992). Note that 
it shows a sign of spiral structure, but without the well- 
defined spiral arms of NGC 3351. Neither does it have the 
expected size and magnitude of a Sbc galaxy at 10 Mpc. 
The radial profile of CGCG 065-091 is shown in Figure 
13b, where the disk and bulge components can be seen. 
Based upon its morphology and properties, this galaxy is 
likely to be in the background of the Leo I group, and this 
fact is corroborated by its velocity. The previously uncat- 
alogued galaxy shown in Figure 14a is a likely non-member 
from the other extreme of parameter space, having a small 
isophotal size as indicated by the detection aperture, and 
having photometric properties (Rt — 19.66, yu = 22.04) 
in Figure 12 outside those expected for LG-type dwarfs. 
The profile in Figure 14b is mostly dominated by seeing, 
due to its small angular size at our isophotal limits, but 
suggests an exponential profile. The profile fit parameters, 
however, are inconsistent with those expected for LG-type 
dwarfs (a m = 1'.'41, nf= 21.95). 

In Figures 15 and 16 we show two objects which are 
each detected by both the traditional and the optimized 
method. In each figure, we show in panel (a) the origi- 
nal MOSAIC image with the detection aperture from the 
traditional method indicated, in panel (b) the masked im- 
age with the detection aperture for the optimized method 
indicated and in panel (c) the convolved image with the 
optimized detection aperture again indicated. A compari- 
son of the aperture sizes in panels (a) and (c) directly con- 
trasts the two detection methods: the detection threshold 
of the traditional method in (a) only detects the bright 
center of the galaxy, whereas the optimized method in (c) 
also detects an extended, low-surface-brightness compo- 
nent. The extremely low surface brightness and extended 
sizes of these objects suggest they are dSphs in Leo I, and 
the object in Figure 15 might be a nucleated dwarf. 

We plot the radial profiles for these two objects in 
panel d of Figures 15 and 16. The profiles extend to 
/Lt(r) ~ 26 R mag/// 2 out to <~ 1.6a fit for Figure 15d, and 
somewhat further in Figure 16d, n(r) ~ 28R mag/// 2 out 
to ~ 3.4a fit . We fit both profiles with an exponential com- 
ponent, adopting ^[j'and R T % . In Figure 15d, where we 
see the suggestion of a nucleus, we fit only the exponential 
component. If the nucleus is real, it should contribute only 
a small amount to the total flux of the object (~ 15% for 
the faintest galaxies) , as has been seen for nucleated dwarfs 
found in the Virgo cluster (Sandage & Binggeli 1984; Im- 
pey et al. 1988). In addition, excluding the nucleus from 
the fit makes the underlying dwarf characteristics more 
directly comparable to dwarf counterparts in the LG (al- 
though, the nature of the nucleus and its relationship to 
the underlying stellar population is still not well known). 
The final parameters of these galaxies support their clas- 
sification as dwarf spheroidal members of Leo I. 

5.2. Supplementary Observations 

Although we have some confidence in the identification 
of Leo-group dwarf galaxy members on the basis of mor- 
phology, total magnitude and central surface brightness 
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(i.e. the location of candidates in Figure 12), we are plan- 
ning and engaged in follow-up observations. To measure 
a completely reliable GLF, we would like to verify mem- 
bership for every Leo member candidate. Particularly for 
the low-surface-brightness candidates, this is an ambitious 
undertaking. One important side benefit of the follow-up 
observations is that when our final list of Leo Group galax- 
ies is assembled, we will have a significant database of the 
structural, photometric, and in some cases, spectral prop- 
erties of the ensemble. We have three on-going programs 
described in the following subsections. 

5.2.1. Multi-color imaging 

Follow-up imaging of our candidates can provide a fuller 
characterization of very low-luminosity galaxies in an en- 
vironment similar to, but distinct from, the LG. In addi- 
tion, the lowest-surface-brightness candidates, which are 
primarily detected through our filtering search technique, 
can be subject to contamination by some false detections 
and high- z galaxies and galaxy clusters. Because of the 
difficulty in flat-fielding the MOSAIC images and the dif- 
ficulty in completely removing chip flaws and the intra- 
CCD dead regions in the final stacked image, some of our 
candidates in the MOSAIC fields are going to be spuri- 
ous. Imaging each candidate with another telescope will 
quickly allow us to cull such detections from our list. 
Our filtering technique also makes us sensitive to back- 
ground LSB disk galaxies (we expect ~ 4/deg 2 , Dalcan- 
ton et al. 1997), cosmologically-dimmed z > 0.4 galaxies 
and z > 0.6 galaxy clusters (< 1/deg 2 ; A. H. Gonzalez 
2000, private communication). These candidates are dif- 
ficult to reach spectroscopically; however, with follow-up 
imaging in B and I, the three colors (we have R from 
the 0.9m+MOSAIC imaging) can flag higher- z galaxy and 
cluster contamination due to their very red colors. In Fig- 
ure 17, adapted from Fukugita, Shimasaku, & Ichikawa 
(1995), we show the distribution of B — R for various mor- 
phological types as they change with redshift (generally 
increasing to the right). The circled values are for z = 0, 
squares are z — 0.2, triangles are z — 0.5 and hexagons are 
z = 0.8. Measured values for LG galaxies today fall along 
the z = locus. Note that all objects with B — R > 1.7, 
according to these models, have z > 0.2. Multiple colors 
can provide a broader baseline for such discrimination, al- 
though such color classification must also be approached 
conservatively, as such models do not account for unusual 
objects, such as the very red LSB galaxies recently discov- 
ered by O'Neil, Bothun, & Schombert (2000). 

5.2.2. Spectral observations 

It is possible to obtain spectra suitable for measuring 
radial velocities for a large fraction of the candidates. We 
have already devoted two Keck nights and two Lick 3m 
nights to this effort with many more required to com- 
plete the follow-up of the Leo Group candidates. Figure 
18 shows an example candidate (2' x 2' postage stamp 
from our MOSAIC image) and its spectrum. In this case 
the spectrum was obtained with the Lick 3m+KAST spec- 
trometer in a 3 x 900s exposure, and the velocity calculated 
via the cross-correlation with a velocity template. The ve- 
locity in this case, 6525 km s _1 , showed this object to be 
background to the group. Figure 19 shows a more chal- 
lenging example, which is an object previously catalogued 



by Ferguson & Sandage (1990) corresponding to LEO #13 
in their catalog list. This spectrum required 1.5 hours in- 
tegration with Keck II+LRIS and cross-correlating with 
a velocity template to extract a velocity. For this object, 
one of the best candidate dSph galaxies, we measure the 
velocity of 886 km s _1 , indicating this is a bona-fide Leo 
Group dwarf. 

5.2.3. Red Giant Branch Tip and Surface brightness 
fluctuations 

At a distance of ~ 10 Mpc (m — M ~ 30) it is pos- 
sible to use two additional photometry-based methods to 
test for group membership for our candidate Leo dwarfs. 
The brightest of the first-ascent red-giant-branch stars, at 
Mi ~ —4, will be measurable at raj ~ 26 mag in Hubble 
Space Telescope images or with large ground-based tele- 
scopes and good seeing. Sakai et al. (1997) used a modest 
amount of HST observing time to demonstrate this for the 
Leo Group giant elliptical NGC 3379. Their observations 
were made in an outer region of the galaxy at a surface 
brightness (fis ~ 27.1, fij ~ 25.2) which is very similar to 
the central surface brightness of our most diffuse candi- 
dates. The tip of the red-giant branch (TRGB) was easily 
detected. 

Measurement of surface brightness fluctuations has 
proven to be a powerful method for estimating distances 
to galaxies to beyond the distance of the Virgo cluster and 
also as a means of inferring the stellar populations in galax- 
ies. NGC 3379 in the Leo I Group has also been used to 
to demonstrate this technique (Tonry & Schneider 1988). 
Bothun ct al. (1991) showed that low-surface-brightness 
dwarf galaxies in the Fornax cluster are particularly well- 
suited for use of this technique and more recently it has 
been applied specifically to low-surface-brightness dwarfs 
(e.g. O'Neil, Bothun, & Impey 1999; Jerjen, Freeman & 
Binggeli 2000). This method has the potential for allow- 
ing membership to be determined for the most challenging 
of the Leo LSB dwarf candidates — those with the lowest 
central surface brightness. Although there is potential for 
learning about the stellar populations in any Leo Group 
members, we will initially be performing a binary test for 
membership. The candidate Leo I galaxies will generally 
either be near 10 Mpc or at a much larger distance in the 
background. The brightest giants will be resolved or not; 
the LSB galaxies will show measureable fluctuation or not. 

In March 2000, we imaged several of the dSph candi- 
dates with the new imaging spectrograph 'ESF on the 
Keck II telescope. With these data we can show the 
promise of the TRGB and SBF techniques for our specific 
application. For three of the galaxies, imaged under excel- 
lent conditions with point sources showing < 0'.'5 FWHM, 
the galaxies just barely resolved into stars in only a lOmin 
exposure in the I band. Figure 20 shows an example of one 
of our 'resolved' candidates (the same object as in Figure 
19). Such excellent seeing conditions are not to be counted 
on for ground-based imaging and we do not expect to be 
able to survey a large number of dwarf candidates with 
this technique from the ground. With HST however, the 
TRGB stars in Leo Group LSB dwarfs can be resolved in 
approximately 20 min observing time. 

It is possible to measure SBF for Leo I dwarfs with 
ground-based imaging. Figure 21 shows the fluctuation 
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power spectrum for the galaxy from the previous exam- 
ple. In Figure 21, we compare the power spectrum from a 
section of the galaxy with no discernable brightness gradi- 
ent (filled boxes) to the power spectrum of a star (stars) 
and to a blank section of the sky (triangles). Note that at 
large k, you see the constant signature of the white noise 
spectrum, whereas at small k, you see the PSF-like signa- 
ture of both the star and the underlying unresolved stars 
in the galaxy image. This contrasts with the spectrum of 
the sky, where there is only a very weak PSF contribu- 
tion from unresolved background galaxies. We have also 
estimated the amplitude of fluctuations via the method 
described in Bothun et al. (1991). The variance in counts 
was measured in multiple 100-pixel boxes placed in 'sky' 
regions of the frame and in six positions within the galaxy 
after point sources were removed and a median-smoothed 
version of the image was subtracted out to remove the 
large-scale gradients in the galaxy profile. The excess fluc- 
tuation signal above that expected in the sky plus galaxy is 
~ 9% of the galaxy signal. We are unambigously measur- 
ing the pixel-to-pixel variations due to the small number 
of average-luminosity giants over the face of the galaxy. 
Rather than attempting to derive an absolute calibration 
for distance based on our SBF measurements in the Leo I 
candidates, we will be able to make a differential distance 
estimate comparing to SBF measured at different surface 
brightness levels in NGC 3379 (and possibly correcting for 
stellar population differences based on color). 

5.3. Preliminary GLF 

With our preliminary membership assignment, we can 
begin to measure the GLF. In Figure 22 we present 
a histogram of absolute i?-band magnitudes (assuming 
in — M = 30 based on cepheid measurements in NGC 
3351; Graham, et al. 1996) for objects found in Field 7 that 
are likely members. While this is not a robust measure- 
ment of the galaxy luminosity density for the group (we 
reserve this analysis to a later paper with the full group 
sample), we note that even with a preliminary member- 
ship classification, we probe fainter magnitudes than have 
been studied in the Leo I group to date. The nominal 
100% completeness limit of Ferguson & Sandage (1990) 
is Ms,iim = —14.2 (Ferguson & Sandage values adjusted 
for m — M = 30), which for B — R = 1.3 corresponds to 
MR t u m = —15.5. Their photographic survey, which covers 
a comparable area on the sky, used an isophotal diameter 
limit of 17" at (ib = 27 mag/// 2 , and beyond this limit 
they applied an incompleteness correction based upon the 
expected range of the fi e — Bt relation as measured in For- 
nax (Ferguson & Sandage 1988). They find this correction 
to be small out to M B - -11.2 (M R ~ -12.5), although 
their completeness description docs not account for objects 
of lower surface brightness for a given magnitude than pre- 
dicted by the [i e — Bt relation like those recently discov- 
ered in Fornax (Kambas et al. 2000). This would suggest 
a larger incompleteness fraction at lower magnitudes, as 
well as perhaps an even steeper faint-end slope than the 
a = —1.4 predicted in Ferguson & Sandage (1991). One 
of the real advantages of our simulation method is that we 
explore these possible surface brightness selection effects. 
We see from our simulations, and from our candidate lists 
that, even at faint surface brightnesses, we can detect ob- 



jects to M R > -11. 

6. DISCUSSION 

Even without full membership determination, we have 
demonstrated the success of our method in discovering 
dwarf spheroidal galaxies in the Leo I group. With the 
advent of MOSAIC imaging, we combine a very wide-field 
survey with the advantages of CCD photometry and can 
now explore nearby galaxy groups to the same sensitivities 
as were previously available only for smaller, and so more 
distant, galaxy groups and clusters. At the 70% complete- 
ness limit for this field, we extend our survey at least 2 
magnitudes beyond the photographic work done by Fer- 
guson & Sandage (1991) for /j,q < 24.5, discovering dwarf 
sphcroidals like that shown in Figure 16. We also find 
an example (Figure 15) of the large, LSB dwarfs first dis- 
covered in Virgo (Impey et al. 1988) — objects which are 
apparently absent from the Local Group. This object cor- 
responds to LEO #5 in the Ferguson & Sandage (1990) 
catalog. 

One possible limitation of our sample selection is that 
our methods are optimized to search for the typical objects 
we would expect in the Local Group. This can lead to some 
prejudices. For instance, very compact, faint objects we 
classify as background, assuming the faintest galaxies will 
be diffuse dwarf spheroidals. Overall, this is justified since 
the only type of compact dwarf seen the local Universe is 
M32-like, and we detect no objects with such high central 
surface brightnesses that SExtractor classifies as galaxies 
(distant objects with these bright are so compact they 
are classified as stars by SExtractor). Furthermore, all 
similarly compact objects detected in the Leo I group by 
Ferguson & Sandage (1990), which were classified as pos- 
sible M32-like candidates but also as non-members, were 
later spectroscopically confirmed as background galaxies 
(Ziegler & Bender 1998). However, unusually compact 
objects not seen in the LG, or BCDs, might be lost from 
our sample in this way. Other objects not seen in the 
LG are the large, LSB dwarfs seen in Virgo (Impey et 
al. 1988), Fornax (Bothun et al. 1991) and M81 (Cald- 
well et al. 1998). In Fornax a population of extreme LSB 
dwarfs was detected that could rival the brighter galaxies 
in integrated luminosity and mass (Kambas et al. 2000). 
Although we don't preferentially search for this type of 
galaxy, we find from our simulations that the sensitivity 
of our search method peaks for bright (i.e. large) LSB 
dwarf galaxies like these. So, if such objects exist in large 
numbers in Leo I (we find only one in Field 7), we will 
detect them. 

Our understanding of the LG also leads us to expect 
dwarfs to cluster as satellites around brighter, companion 
galaxies. This correlation is seen in some other low den- 
sity environments (Vader & Sandage 1991; Loveday 1997), 
although there are hints that the dwarf-giant correlation 
is lowest in poor group environments, where dwarfs might 
even be 'free-floating' (Vader & Chaboyer 1992). To ex- 
plore this, we have tried to recover as much area as possible 
surrounding the large galaxies in our sample, but for late- 
type, L* galaxies like NGC 3351 in Field7, there would 
be significant confusion of diffuse dwarfs with any residual 
structure from the profile subtraction, such as HII regions 
and spiral arms. Thus, we ignore any detection within a 
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projected distance of <~ 4' (11.5 kpc) around NGC 3351. 
Despite this limitation, however, we would still expect to 
detect a significant fraction of a satellite dwarf population. 
The characteristic distance of a Milky Way dwarf satel- 
lite is at a projected distance of 141 kpc from the Galaxy 
(Karachentsev 1996), which would correspond to ~ 48' at 
10 Mpc, and the corresponding distance for M31 satellites 
is 85 kpc (Karachentsev 1996), or ~ 30' at 10 Mpc. As 
we do not see the Milky Way satellites in projection, it 
is difficult to quantify how many would fall too close to 
the Galaxy to be detected. In the case of M31 satellites, 
two objects would fall within 4' of the main galaxy, but 
these arc M32 and NGC 205 which cause little concern for 
confusion. 

The primary limiting factor for these data is the nature 
of the engineering-grade CCDs used in the MOSAIC cam- 
era for the majority of our Leo I images. Even modest 
brightness enhancements due to lack of flatness of the im- 
ages can be enhanced through our filtering process, forcing 
us to raise our detection threshold to exclude possible spu- 
rious detections. With flatter images, we can lower our 
threshold and push our limits to fainter surface bright- 
nesses and magnitudes. MOSAIC images of our Fields 8 
and 9 were taken with the newer science-grade CCDs last 
Spring, and will provide this comparison. Our current lim- 
itations are well quantified, however, with our simulations. 

The galaxy simulations not only highlight selection ef- 
fects in detection, but also our ability to measure photo- 
metric parameters for these candidate galaxies. As we 
have shown in §4.3, parameter recovery can be rather 
difficult for the lowest-luminosity and lowest-surface- 
brightness objects. Part of this difficulty is in using an au- 
tomated process, and for all real objects detected through 



the optimized method, we opt to fit them by hand. Full 
understanding of these challenges to our photometry pro- 
vides a robust and direct measurement of our errors. 

The discovery of dwarf spheroidal galaxies outside the 
Local Group is a burgeoning pursuit that has begun to 
yield tantalizing results in a number of environments (e.g. 
Virgo and Dorado clusters; Phillipps et al. 1998; Carrasco 
et al. 2001). With theoretical models predicting a steep 
faint-end slope to the GLF (Cole et al. 2000) and detec- 
tions of possible bright dSph progenitors seen at z > 0.2 
(Guzman et al. 1998), it is critical to establish a census 
of these low-luminosity galaxies at z = 0. As we clearly 
establish the selection parameters influencing our sample, 
we can construct the only complete sample of dwarf galax- 
ies to these luminosities and surface brightnesses outside 
of the Local Group. Furthermore, with such a sample we 
may explore the intrinsic nature of these galaxies as an 
ensemble, such as their spatial distribution, abundances, 
and stellar populations. 
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Fig. 1. — Upper: The typical ranges of i?-band central surface brightness vs. total absolute magnitude occupied by local spirals, ellipticals, 
globular clusters, irregulars and dwarfs, adopted from Ferguson & Bingelli (1994) Figure 3. Plotted points arc the members of the Local 
Group, taken from Mateo (1998, assuming B — R = 1.3 for dwarfs without a R magnitude) and in the case of And V, And VI and And VII 
taken from Caldwell (1999). For the Milky Way and M31, the central surface brightnesses arc from Freeman's Law (Freeman 1970). Local 
Group galaxies are coded by morphological type: five-point stars are Spirals, open triangles are Irrs, filled triangles are dlrrs, circled filled 
triangles are transition dSph/dlrr and filled circles are dSphs. The parameters of M32 and Malin 1 are plotted individually to demonstrate 
their extremity. Lower: A blow-up of the region of parameter space that our Leo I survey will probe. Local Group galaxies are now plotted 
with central surface brightness vs. apparent J?-magnitude as they would be seen at the distance of Leo I (10 Mpc). Symbols arc the same as 
in the upper panel. The solid lines represent lines of constant angular size at an isophotal limit of fi = 26.7mag/// , assuming an exponential 
profile (Allen & Shu 1979). Dashed lines arc lines of constant exponential scale length, for a purely exponential profile. 

Fig. 2. — Our coverage of the Leo I group, overlaid on an image of Leo I from the DSS. We show our nine 1° X 1° fields imaged in R with 
the KPNO 0.9m+MOSAIC. For reference, NGC 3379 can be seen near the center of Field 4. We use Field 7 in an example of our search 
technique and first results. 

Fig. 3. — A 15' X 15' sub-image of Field 7 with artificial galaxies added that mimic the characteristics of Local Group dwarfs if seen at a 
distance of 10 Mpc. Note that Leo II, Tucana, And V, and Carina would all be difficult to detect by eye. 

Fig. 4. — The results of Step 1 of our detection method on the image from Figure 3. A SExtractor aperture is drawn around every object 
detected that has been classified as a galaxy (see text). Most artificial dwarfs are detected, except for galaxies with the characteristics of And 
II, Tucana, And V and Carina, which have central surface brightnesses fainter than the detection threshold. 

Fig. 5. — The results of Step 2 of our detection method, where we show the masked and convolved version of the image from Figure 3. All 
SExtractor detections arc indicated, where obvious masking remnants from stars or bright galaxies are marked with crosses and candidate 
detections are marked with SExtractor apertures. Note that in this method, And II, Tucana, And V and Carina are now robustly detected; 
however, we lose Leo II due to its center being masked. This demonstrates the necessity of using both methods together. 

Fig. 6. — Total apparent magnitude vs. central surface brightness for the detections from Figures 4 & 5. Detections from the traditional 
method arc shown as crosses and detections from the optimized method as stars. Local Group dwarfs as seen at 10 Mpc arc plotted as filled 
circles. The inputted artificial galaxies from Figures 4 & 5 are labeled and connected by a line to their measured parameters from both 
methods. The solid and dashed lines are the same as in Figure 1. 

Fig. 7. — a: Detection efficiency contours from our simulations of recovering input artificial dwarfs in Field 7. These contours are for objects 
recovered only by the traditional method. The greyscale shows recovery fractions of: 90% (darkest), 70%, 50%, 30%, 10%(lightest). Local 
Group galaxies if seen at 10 Mpc arc again plotted as filled circles, and the solid and dashed lines follow Figure 1. b: Detection efficiency 
contours from the simulations, showing the efficiency of detecting objects that are only found with the optimized method. Greyscale, symbols 
and lines are the same as in Figure 7a. c: Detection efficiency from simulations for the combined methods. Greyscale, symbols and lines are 
the same as in Figures 7a & 7b. 

Fig. 8. — Parameter recovery as measured by the simulations on Field 7, showing the difference in output and input central surface 
brightness (/Jq u * and n'J 1 , resp.) as a function of input total magnitude, Rip. The output central surface brightnesses are extrapolations from 
exponential fits to the radial profiles. The left panels are all the objects detected through the traditional method, while the right panels are 
the object only detected through the optimized method. Objects are also divided by input central surface brightness in the vertical direction. 

Fig. 9. — Recovery of total magnitude as measured by simulations, where the output total magnitudes (ffip 1 ') for both detection methods 
arc MAG_BEST from SExtractor. The difference in output and input magnitude is plotted as a function of input total magnitude (Rip). The 
divisions by panel are the same as in Figure 8. 

Fig. 10. — Recovery of total magnitude, as in Figure 9, except the output total magnitudes plotted for both methods are now calculated 
from the exponential profile fits. Panel divisions are the same as in Figure 8. 

Fig. 11. — Recovery of exponential scale length (a) from simulations, where output exponential scale length (a out ) is measured from the 
exponential profile fit of recovered objects. Again we plot the difference in output and input a as a function of input total magnitude. Panel 
divisions are the same as in Figure 8. 

Fig. 12. — Same as Figure 6, except here we plot all detections from Field 7. Objects detected by both methods arc connected with a line. 
Example objects discussed further in §5 are circled. 

Fig. 13. — Detected background spiral, CGCG 065-091. a: 2' X 2' postage stamp of the MOSAIC ij-band image, with detection aperture 
from SExtractor indicated, b: Radial profile with errorbars as measured by ELLIPSE in IRAF. Fitted exponential disk and Ri-Rt overlaid. 

Fig. 14. — Detected distant background objects — probably a distant spiral, a: 2' X 2' MOSAIC ij-band image, as in Figure 13. b: Radial 
profile of object, as in Figure 13. 

Fig. 15. — [figurel5.gif] Candidate nucleated dwarf, a: 2' X 2' MOSAIC R-band image, showing the traditional method detection aperture, 
b: The same image, masked of bright objects, with the optimized method aperture overlaid, c: The same image as in b, convolved, with the 
optimized aperture overlaid, d: The measured radial light profile with exponential fit drawn as a line. Fitted exponential scale length, ct flt , 
and half-light radius, r^g, are labeled. 

Fig. 16. — Newly discovered dwarf spheroidal. Panels same as in Figure 15. 

Fig. 17. — A diagram of model B — R values from Fukugita ct al. (1995), which shows the distribution of color for each morphological type 
as it changes with redshift. The circled values are for z = 0, squares are z = 0.2, triangles are z = 0.5 and hexagons are z = 0.8. Measured 
values for Local Group galaxies today correspond to values along the z = locus. Note that all objects with B — R > 1.7 have z > 0.2. 

Fig. 18. — Candidate galaxy found to be in the background of Leo I. a: 2' X 2' image as in Figure 13a. b: 3x900s spectrum taken with 
Lick 3m+KAST. Velocity measured is 6525 km s _1 , which is far in the background of the group. 

Fig. 19. — Candidate dwarf confirmed to be in Leo I. a: 2' X 2' image from MOSAIC. Labeled parameters are from profile fit. b: 1.5 hour 
spectrum with Keck II+LRIS. Measured velocity is 886 km s -1 , placing this dwarf well within the Leo I group. 

Fig. 20. — /-band image of Leo I candidate dwarf shown in previous figure. This image was taken on Keck II+ESI with extraordinary 
seeing (< 0'.'5 FWHM), and quickly resolved into the brightest stars. Removing the point sources from this image, we can measure the 
surface- brightness-fluctuation signal. 
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Fig. 21. — The SBF power spectrum of the object in Figure 20. The filled boxes are the power spectrum measured from a patch of the 
galaxy, the stars arc the spectrum measured from a bright star in the image and the triangles are the spectrum we measure from a randomly 
chosen empty region of sky. The galaxy spectrum (boxes) shows we are measuring a PSF-like signal from unresolved stars in the galaxy which 
traces the PSF spectrum from the star (stars) and far exceeds the very weak PSF-signal from unresolved background galaxies in the empty 
sky image (triangles). 

Fig. 22. — A histogram of detected objects in Field 7, applying a preliminary morphological membership classification. While this is not 
yet a robust measurement of the Leo I luminosity function, it demonstrates the depth of our survey. 
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